Intergalactic cold dust in the NGC 4631 group 
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Abstract 

We have detected extraplanar cold dust at distances out to > 10 kpc, situated in the halo of the interacting 
galaxy NGC 4631. The dust emission disk is much thinner than the warped Hi disk and new structures 
emerge. In particular, a giant arc has been found that is linked to anomalies in the kinematical structure 
of the atomic gas. Most of the extraplanar dust is closely associated with Hi spurs that have been found 
earlier [Weliachew, L., Sancisi, R. & Guelin, M. (1978) Astron. Astrophys. 65, 37-45; Rand, R.J. (1994) 
CO ■ Astron. Astrophys. 285, 833-856]. These spurs obviously are traces of the interaction [Combes, F. (1978) 

Astron. Astrophys. 65, 47-55]. The dust emission within the plane reaches the border of the optical disk. 

The activity of the disk of NGC 4631 is moderately enhanced by the interaction, but no gas moving 
in the z-direction could be found [Rand, R.J., Kulkarni, S.R. & Hester, J.J. (1992) Astrophys. J. 396, 
97-103; Golla, G., Dettmar, R.-J. & Domgorgen, H. (1996) Astron. Astrophys. 313, 439-447]. Hence, it 
seems unlikely that strong winds have deposited the high-z dust. Instead, the coincidence with the Hi 
features suggests that we see a track left behind by the interaction. In addition, the Hi shows a supershell 
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formed by an impact [Rand, R.J. & Stone, J.M. (1996) Astron. J. Ill, 190-196] in the zone where the 
dust trail crosses the disk. This region is also characterized by disturbances in the distribution of the 
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Ha light. The masses associated with the dust can be estimated only very roughly on the basis of the 
existing data; they are of the order of a few 10 9 M of gas. 



^ ■ 1 Introduction 



Cold dust has come into focus only recently because it had to await the development of sensitive 
millimeter /sub- millimeter bolometer arrays to be detectable unambiguously. The IRAS survey could 
provide only hints at its existence because it was blinded by the strong emission from the small percent- 
age of warmer dust that is radiating far more brightly The large amount of cold dust (Td < 25 K) can be 
detected only at (sub)mm wavelengths where the radiation of the warmer components has vanished. To 
give an example: The peak brightness of a blackbody at 30 K is 30 x higher than that of a 15-K object, 
all other parameters being equal. On the other hand, the radiation of a blackbody at 15 K peaks at 
~ 200 fim and remains more than one order of magnitude brighter at A 1.2 mm than that of a blackbody 
at 30 K with the same peak brightness. Now, the emissivity of interstellar dust is roughly proportional to 
T® - for a blackbody B(T) = a T 4 - so even a very large amount of cold dust emits only weakly. Because 
of this T 6 dependence of the emission, a very large energy input is needed to heat dust, and the majority 
of it remains at lower temperatures. This cold component thus is an important tracer, and, indeed, it 
may represent > 90% of the interstellar dust (cf. refs. 1 and 2). 

In itself, the contribution of the dust to the total mass of a galaxy is < 1% of the gas mass, but there 
are indications that the dust-to-gas ratio is relatively constant, independent of the type of gas (atomic or 
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molecular). Indeed, the dust grains are believed to play a crucial role in the formation of molecules. The 
study of their properties thus also offers an independent means of studying the molecular gas content of 
galaxies. This is important because the standard practice of observing the CO molecule and deriving, 
thereby, the properties of the H2 has substantial uncertainties particularly concerning the derived masses. 
On the other hand, investigating the cold dust is technically difficult and cannot provide any information 
about the kinematics because it is based on broadband continuum observations. 

In several runs, the IRAM (Institute for Radio Astronomy in the Millimeter domain) 30-m telescope 
equipped with MPIfR (Max-Planck-Institute for Radioastronomy) bolometer arrays has been used to 
map nearby galaxies in the A 1.2-mm continuum emission. The first maps led to the impression that 
it is well correlated with the CO emission, and drops off similarly steeply with increasing distance from 
the center. This behavior was shown, for example, for galaxies NGC891 (3), M51 (4), and NGC4631 
(5). It soon became evident, however, that this is not generally the case. The Sb galaxy NGC4565 is 
significantly more extended in the emission of A 1.2-mm continuum than in that of the CO line (6). The 
cold dust is even detected in the warped outermost rim of the disk. As an intermediate case, NGC 5907 
also shows an extended dust disk (7). The sensitivity needed to detect this extended emission has been 
achieved only recently, however, and the sample is still small. So it is not yet clear what determines the 
extent of the cold dust - the profile of the cold dust along the major axis in NGC 891 remains very close 
to the rapidly vanishing CO, even when studied with much higher sensitivity than previously published 
(R. Zylka, personal communication). 

2 The observations and the object: NGC 4631 

2.1 Observational details 

All recent maps were obtained with bolometer arrays consisting of 19 elements whose sensitivity is in 
practice about a factor of 2 better than that of the 7-element detector used before. The 19-element 
bolometer array has a bandwidth of ~ 80 GHz centered at ~ 230 GHz. The individual elements are 
arranged in a closely packed hexagonal pattern. The beam size at the 30-m telescope is 11" and the 
spacing between the beams 20". The observations were made in March 1997 during a period of stable 
weather with zenith opacities typically < 0.2. We monitored the sky opacity before and after each 
subimage and mapped Mars every night to determine the absolute flux scale. To obtain a map, the 
object is scanned in azimuthal direction including parts of blank sky on both sides to define a proper 
zero level. In addition, the sub-reflector of the telescope is oscillating at a frequency of 2 Hz which makes 
the beam switch between two positions separated by 45" in the orientation of the scanning. This yields 
an "on-off" measurement that cancels atmospheric variations at short time scales. The whole area of 
NGC 4631 was covered with a mosaic of 19 individual fields. Their distribution is rather uniform along 
the whole disk and, hence, the sensitivity drops only at the outer edges of the final map. Furthermore, we 
kept the scanning orientation close to the minor axis of the galaxy by carefully chosing the hour angles of 
the individual observations. This minimizes spurious contributions and the noise in the map. The field 
presented here is about 15' x 8' after the cutoff of the edges with lower sensitivity. In the central part, 
the noise is ~ 2 mJy/beam for the data smoothed to an angular resolution of 20" and rises to ~ 3.5 
mJy/beam at the edges. 

2.2 NGC 4631 

It is obviously best to choose edge-on galaxies for studies of weak phenomena, because the lines of sight 
are long - remember that the dust emission in the mm regime is optically thin and thus the whole disk 
contributes to the detectable flux. The 19-element bolometer made it possible to map large areas, so we 
decided to re-observe NGC 4631. This moderately active galaxy has long been a favorite candidate for an 
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interacting system. It is relatively nearby, at 7.5 megaparsecs (Mpc) (ref. 8; 1' corresponds to ~ 2kpc), 
and two obvious companions are close by. The dwarf elliptical NGC4627 is situated 3' northwest of 
the nucleus and 30' to the southeast the distorted spiral NGC 4656 can be found. The whole group has 
been extensively studied in Hi (8,9) in order to understand the traces of the interaction (cf. Fig. ^). 
According to a modeling of the encounter (10) the prominent streamers of atomic gas can be explained 
as being pulled out of the members of this group during the interaction. Presumably also as a result of 
the interaction, the disk of NGC 4631 has a disturbed appearance in the optical continuum and Ha line 
emission. At two positions in the disk highly energetic "supershells" have been found (11) of which one 
is described as being caused by the impact of a high velocity object (12). 

Almost simultaneously with the early high-resolution Hi observations a large radio halo was detected 

(13) . It is of nonthermal origin and one of the most prominent radio halos known. Detailed investigations 

(14) have subsequently shown that the magnetic field lines of the galactic disk open into the halo - in 
stark contrast to most spiral galaxies, where the field is more or less confined in the disk (15). This 
magnetic field structure allows electrons, cosmic ray particles and hot gas to escape from the active disk 
into the halo. ROSAT (Roentgen Satellite) detections (16,17) show a large X-ray envelope that is a 
natural consequence of this configuration. 

The investigation of the molecular gas indicates rather normal conditions, however. The central region 
of ~ 2' x V was completely mapped in the (1-0) and (2-1) transitions of CO (18). In addition, a major 
axis strip of 7' length was obtained with a higher sensitivity of ~ 20 mK. Some additional spectra at 
other locations did not reveal significant emission. 

3 The distribution of the cold dust 
3.1 The dusty disk 

At a first glance, the A 1.2-mm emission of NGC 4631 (Fig.0) is characterized by a narrow, extended 
disk, with a double-peaked central region. About three times weaker than the brightest peaks, at a level 
of ~ 30mJy/beam, the disk is stretched out to a distance of ~ 13kpc on either side of the nucleus, 
gently decreasing in the west, and at a relatively constant level for some lOkpc in the east - see Fig.||. 
This distribution is similar to that found in NGC 4565 (6): The correlation between dust and CO is 
restricted to the nuclear region, whereas the dust in the outer parts of the disk seems to follow the Hi. 
The CO emission drops in a similarly steep manner as the centimetric radio continuum: At a radius of 
2.5' (= 5.5 kpc) it has decreased to a tenth of the peak value, and a bit further out, no CO emission 
could be detected at all. At this radius, the dust emission is still at about one-third of the peak level, 
however, and it actually stretches out at least to the edge of the optical disk. In fact, the limits are set 
by the border of the map, not by the vanishing emission. The old bolometer map (5) is limited to the 
innermost part due to its restricted coverage and sensitivity. 

The narrow main emission ridge is somewhat surprising given the fact that NGC 4631 does not even 
show a dust lane. Its thickness is similar to that of the undisturbed edge-on galaxies NGC 4565 (6) or 
NGC 891 (3). This suggests that the dust is concentrated in the midplane of a galaxy in any case (cf. 
ref. 19), and the absence of an optical dust lane may just reflect a high dumpiness with large "holes". 

A closer inspection shows some radial variations in comparison with the Ha map (Fig. |^). Near 
the center, the optically bright regions surround the dust emission peaks. Here, the bulk of the optical 
emission is produced in the central region and later absorbed on its way through the disk. In contrast, 
in the outer parts of the disk several dust maxima coincide with optically bright spots. It is unlikely 
that all of these places are just accidental line-of-sight coincidences. The source of the radiation could be 
young stars in their dusty birthplaces (e.g. in a spiral arm tangent). In any case, the material along the 
light path should be rather transparent. Such large variations of the opacity are not surprising for spiral 
galaxies, however (see ref. 20 for a compilation). 
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3.2 The extraplanar dust 

Completely new is the detection of significant dust emission out to z-distances of at least lOkpc. The 
distribution of this intergalactic dust strongly suggests that it has been brought there by the same 
mechanism that formed the four Hi streamers (8,9). All three Hi spurs that are touched by our map are 
connected with the thin optical disk by corresponding dust features (see Fig. ||). The disk formed by the 
atomic gas is much thicker than the optical or dust emission disk, however, so that most of the A 1.2-]mm 
emission still lies within their boundaries. On the other hand, CO emission could only be detected out 
to z ~ 1 kpc. So either there is hidden molecular gas in those outer regions or the dust is associated with 
atomic gas. 

Although these streamers follow the already known Hi features rather closely, north of the center of 
NGC 4631 a structure has been unveiled that was invisible in the thick atomic gas disk: A giant arc spans 
over the central region with its footpoints about 4 kpc east and west of the central region. The eastern 
footpoint is situated opposite the onset of the southern streamer 2, and the western part of the arc is 
blending into spur 4. The thickness of the disk emission seems to be reduced in the central region, but 
this might be an artifact of the data-reducing technique. The question arises regarding whether there is 
a common origin for these extraplanar structures. 

4 Origin of the extraplanar dust 

Unfortunately, we do not have any velocity information from the A 1.2-mm continuum observations, and 
hence have to rely on indirect arguments for the determination of the history of this dust distribution. 
To some degree, it seems reasonable to assume that the Hi velocities are a good indicator also for the 
dust kinematics. Indeed, there are at least two clear coincidences between Hi velocity anomalies and the 
arc. The velocity gradient along the major axis is much steeper in the northern part of the disk than 
in the southern. A possible interpretation suggests that the disk is somewhat inclined and warped along 
the line of sight, so we would be looking at different portions of the galaxy (8) - north of the center at 
the inner disk with a steeper gradient and south of it at the near part of the outer disk. 

If we compare the velocity field with the dust emission (see Fig. ||), a different explanation arises. In 
addition to the steeper gradient, the northern disk shows anomalous velocity components at two places 
- and they coincide perfectly with the footpoints of the arc. These anomalous components are visible as 
regions with almost closed isovelocity contours in Fig. |s| (marked as Al and A2). Although we are not 
able to determine the precise location of the extragalactic dust with respect to the disk, it seems clear 
that there is a local interaction. The additional velocity components point towards us in the west and 
away from us in the east. 

In addition there are two structural indications, both of them suggesting that material has followed 
the arc in an counterclockwise direction and hit the disk in the east: (i) The appearance of the disk in the 
light of the Ha line is much more disturbed at the eastern footpoint of the arc (cf. Fig. [| In broadband 
red light the eastern part seems truncated. The disturbed region lies between the arc and the streamer 
2, whereas the western side seems unaffected by the arc or streamers 3 and 4. (ii) In the Hi emission 
of NGC 4631 two supershells have been found near the midplane and subsequently have been modeled. 
Shell 2 (west of the nucleus at Right Ascension 12 ,l 39 m 30 s ) seems to be an expanding bubble, whereas 
shell 1 (at 12' l 39 m 56 s ) is most probably caused by an impact of a cloud with a mass of ~ 10 7 M© at a 
speed of 200 km/s coming from the north (12). The geometry of the collision is very well constrained by 
its kinematical signature in the Hi data and indicates that he material not only came from the north but, 
moreover, must have had a velocity component away from us. The location of this shell is between arc 
and spur 2. This fits to the other indications of the trajectory of the arc (the shells and the anomalous 
components mentioned above are distinct features). 

If we try to describe all extraplanar dust features as a single trail, its path could be as follows: starting 
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from spur 3, it runs through the western disk, follows it to the western footpoint of the arc, sweeps along 
the arc, penetrates the eastern disk and leaves it via spur 2. In the way, some material is forced away: 
for example, to follow streamer 4. Of course, such a concatenation is purely speculative and, at this 
point, merely meant to summarize the structure of the dust distribution. In principle, the concept of a 
continuous trail is a natural explanation within the framework of an interaction, however. In any case, 
we have to answer the question regarding which of the three involved galaxies has left the dust traces. 

The interaction has been modeled rather early on the basis of the first high-resolution Hi data (10). 
This description suggests that the bridge 1 and the streamer 4 had been parts of the disk of NGC4631. 
The two perpendicular features 2 and 3 arc made from material pulled out of the now dwarf elliptical 
NGC 4627 which might have been of a different type before. However, the Hi emission of structures 2 
and 3 is only weakly connected with the disk of NGC 4631 and not at all with the dwarf. It has to be 
stressed that this model is by no means unique - three involved galaxies open a vast parameter space for 
an interaction scenario. 

Today, the location of part of the Hi gas can be determined even in the third dimension with the help 
of X-ray data: The soft band is very sensitive to absorption by atomic hydrogen, and it can be safely 
concluded from the distributions that the southern part of the disk is at the near side (16,17). Looking a 
bit more in detail, the onset of spur 4 seems to be located in front of the X-ray halo (cf. Figure 6 of ref. 
17) but, higher up, the X-ray emission becomes stronger again, suggesting that this streamer is pointing 
away from us. Such additional information will help to constrain the parameters of the interaction much 
better than it was possible before. 

5 Discussion and outlook 
5.1 Origin of the mm emission 

NGC 4631 is characterized not only by the interaction with two other galaxies but also by a unusually large 
radio halo. In particular the enhanced star formation might be responsible for the uncommon magnetic 
field configuration and thus for the radio and X-ray halo - could this possibly produce extraordinary 
continuum emission at A 1.2 mm as well? 

In addition to radiation from dust, there are three candidates for the source of radiation: free- free 
emission in ionized clouds, nonthermal emission and line radiation within the bandpass of the detector. 
The free-free emission should be correlated with Ha - so a contribution is only expected in the disk since 
the sparse high-z Ha emission (21) is not correlated with dust features. Within the disk, even in the 
brightest spots the emission measure reaches only a relatively low value of 1000pc-cm~ 6 (22). So the 
total flux due to free- free radiation is of the order of lOmJy only. 

The synchrotron emission is very extended around NGC 4631 at cm wavelengths (14). From the fluxes 
and the determined spectral index (23) we can estimate its contribution in our mm band. In the center, 
the spectral index is rather flat at about a ~ —0.65. The peak flux here is 60 mJy/84" beam at 10.55 GHz, 
this translates into 0.3 mJy/20"beam at 230 GHz. Outside the disk the spectral index is even steeper, 
so this contribution is smaller than a tenth of the noise level. The only significant addition thus comes 
from the molecular line emission: The eastern peak reaches 70 K-kms -1 of 12 CO(2-l) emission. Within 
our 20" beam, this contributes to a flux density of ~ 20mJy/beam (cf. ref. 4), about a fifth of the total 
value. In the outer disk, however, the line emission adds < 3mJy/beam - just the figure of the noise 
level in the map. 

So, the observed mm continuum radiation is pure thermal dust emission almost everywhere, but how 
did the dust reach such enormous z-heights? Could it simply be blown out of the active plane by large- 
scale winds? Outflows are known for many galaxies with strong starbursts like, for example, M 82 (see, 
e.g., refs. 24 and 25). But the disk of NGC 4631 is forming stars at an only moderately enhanced rate 
(18). Ha kinematic data as a more direct tool do not show signs of a gas outflow from the disk (22) - in 
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contrast to M82 which has a similar orientation (24). Peculiar velocities have been found, but they are 
more likely explained by the direct influence of the interaction. 

5.2 Temperatures and Masses 

The dust around the disk of NGC4631 seems to be cold - the spatially resolved IRAS CPC data as an 
indicator for warmer dust show at most a somewhat thicker disk near the center (26). Its width is of 
the order of 90" at lOO^m and 45" at 50/xm wavelength using the published beam sizes of 95" and 80", 
respectively. No significant emission can be seen further out. Unfortunately, the CPC instrument could 
not be calibrated properly, and because of the uncertainties of the order of ±60% we could not derive 
dust temperatures from these data. 

In view of the moderate star forming activity the presence of very warm material does not seem 
very likely. Presumably the temperature of the coldest dust component is rather low even in the disk - 
similarly to other spirals in which it could be measured so far (2,6). We expect, therefore, temperatures 
in the range of 15 to 20 K for the disk, and even lower values for the extraplanar dust. To determine such 
temperatures, observations in the sub-mm range between 1 mm and 100 /im are needed, but, here, the 
opaque atmosphere renders them very difficult. This is why only few data of spiral galaxies have been 
published in this range (e.g. refs. 2 and 27). We will do more observations at AA 450 /im and 850 /im 
to complement the existing data. Together, they should provide the spectral information needed for the 
determination of the temperatures. 

Until then, we have to postpone the proper determination of masses and energies from the dust 
emission. In principle, the dust mass is directly proportional to the detected flux S\. With some 
additional assumptions about the dust and gas properties, we can derive the hydrogen column density. 
In the millimeter/sub- millimeter regime, the observed flux density per beam (of half power beam width 9) 
produced by dust of temperature Td is given by S\ = 1.13 9 2 (1 — e Tx ) B(X, Tj). Here, t\ — <x^ • Ah gives 
the dust absorption cross section per hydrogen atom, and B(X, Tj) is the Planck law for the radiation of 
a black body (see ref. 1 for a derivation, and ref. 6 for a determination of a possibly typical value). Such 
a calculation yields about 2 • 10 9 M© of hydrogen for the gas associated with a dust component of 21.5 K 
in the central region (5). As already stated above, the warmer dust radiates much more efficiently, and, 
hence, a second dust component at 55 K is associated with < 1% of this mass (5). 

These values are poorly constrained, however, and we face large uncertainties in the temperatures of 
the coldest component. If we nevertheless assume 15 K for the extraplanar dust the total gas mass in the 
arc should be of the order of 1.7 • 10 9 M Q , somewhat less than the mass in the central region. A similar 
value is found for its presumed continuation, spur 2, so that the total gas mass in this structure would 
be about half of the atomic hydrogen mass of the whole galaxy (8). 

5.3 Conclusion 

Obviously, the interaction is the dominant event in the history of NGC4631. It has caused a partial 
disruption of the disk, triggered the star formation activity, provoked an upturned magnetic field con- 
figuration, and, last but not least, left trails of atomic gas and dust all around the disk. It is therefore 
crucial to know exactly how the interaction took place. The existing model (10) was obtained by using 
a "trial-and-error" method for the parameters. It is virtually impossible to find the proper solution that 
way; the number of possible parameter sets for an interaction in such a group with three members is 
simply too large. Moreover, the now existing data (including the distribution of the cold dust) give new 
and more detailed information to check the model against. The problem is how to explore the parameter 
space to find the best set. Now there are new, powerful techniques available like so-called 'genetic codes' 
(28), so we hope to obtain a much better view of the interaction, which in turn should help to fix the 
roots of the related phenomena. 
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In any case, the detection of cold dust outside of a galactic disk has unveiled the existence of hitherto 
invisible baryonic matter in the halo region. This will not solve the dark matter problem, even using 
a favorable estimation of the associated mass, but it is a step towards a more "normal" neighborhood 
of galaxies - containing less exotic material, probably simpler to understand, but in any case easier to 
investigate. 
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Figure 1: Map of the A 1.2-mm emission of NGC4631, overlaid on an image taken from the Digital Sky 
Survey. The levels are -6 (dotted), 6, 11, 21, 41, and 81 mJy/beam. Only significant emission is shown 
and the outer parts of the map with higher noise have been cut off. The small object north of the disk 
is the dwarf elliptical galaxy NGC 4627; the other companion, NGC 4656, is situated about half a degree 
away in the south-east. 
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Figure 2: Cuts along the major axis; the CO and the 2.8 cm data are provided by G. Golla (14,18) and 
the Hi curve is extracted from a map supplied by R. Rand (8). The left axis gives the units for the dust 
and CO emission, the right axis for Hi and cm continuum. Similar to the galaxies NGC 4565 (6) and 
5907 (7), the profiles of the CO and the cm continuum drop off more steeply than that of the dust and 
the Hi. Note, however, that such cuts are one-dimensional and miss emission that is close to, but not on 
the major axis. 
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Figure 3: Comparison of the A 1.2-mm emission of NGC4631 with the Ha emission which has been 
smoothed to 6". The dwarf elliptical galaxy NGC4627 is invisible here. Note the large disturbed area 
in the west with its subsequent 'cutoff'. Near the center, the dust and Ha tend to be anti-coincident 
whereas in the outer disk there are clear correspondences. 
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Figure 4: An overview of the distribution of the extraplanar Hi (purple contours) and cold dust (in color); 
the Hi is smoothed to 22" spatial resolution. For comparison, optical isophotes are added in red. The 
Hi spurs are numbered as introduced in ref. 9. The axis labels are given in arc minutes from the center 
of NGC 4631 on the lower and left sides; the other sides gives the projected length scale for an assumed 
distance of 7.5 Mpc. 
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Figure 5: A comparison between the cold dust emission (at 20" resolution) and the velocity field of the 
atomic gas (at 12" x 22"). Velocity contours are plotted every 20km/s (cf. also the right-hand scale). 
Note the velocity anomalies at the foot points of the arc (arrows). The insert shows position- velocity 
curves parallel to the major axis; the box is labelled in velocity vs. minutes of arc. The red line goes 
through Al, the green line through A2 as indicated by short lines in the main plot. In the southern disk 
the gradient is flatter (purple line). The stars mark the position of the two Hi supershells. 
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